Nonrelativistic ionization energy for the helium ground state. 



(N 
O 

o 

o 
O 



Oh! 



o 

-I— > 
03 

c/5 
O 

>V 
^ ■ 

(N 
> 

o 
o 

O 

o : 
• i— i 

^ : 

Ph. 



Vladimir I. Korobov 
Joint Institute for Nuclear Research 
Ul 980, Dubna, Russiaf\ 

The helium ground state nonrelativistic energy with 24 significant digits is presented. The calcu- 
lations are based on variational expansion with randomly chosen exponents. This data can be used 
as a benchmark for other approaches for many electron and/or three-body systems. 

PACS numbers: 31.15.Pf, 31.15.Ar 

From the early days of the quantum mechanics the ground state ionization energy of the helium atom was a 
benchmark for approximate methods of solving nonrelativistic Schrodinger equation for a few-body system. One of 
the earliest variational calculations has been performed by Hylleraas M in 1929 and it yields an about 5 significant 
digits for the nonrelativistic ionization potential. In 1957, Kinoshita p[ presented a 7 digit number obtained with a 
39 parameter function, which along with higher order corrections including the Lamb shift calculations confirmed a 
very good agreement with the best experimental value. Since that time with the development of computer power the 
accuracy grows very rapidly. We would like to mention here the two most recent calculations. The first p| is aimed 
to elaborate an efficient variational method for the many electron atoms. The second is to find an effective and 
economical way for studying the helium and helium-like two electron atoms. 

In this short report we want to present a new very accurate value for the nonrelativistic energy of the helium 
ground state. In our calculations we strictly follows a method described in ||. The two modifications are to be 
stated. First, a sextuple precision arithmetics (about 48 decimal digits) implemented as a Fortran-90 module has 
been used instead of a multiprecision package written by Bailey. The use of this module gives an about 5-fold 
improvement in computational time and allows for to increase significantly a length of a basis set. The module is 
based on a representation of a sextuple precision number by a set of 3 double precision numbers. It is assumed that 
an exact sum of these double precision numbers is some sextuple precision number Second, we have taken a 
multilayer variational wave function with 5 independent sets of variational parameters (instead of 4 as in ||) , which 
consecutively approximates one after another smaller and smaller distances of electrons with respect to a nucleus. 

A variational wave function is expanded in a form Q 
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where a* , Pi and 7, are complex parameters generated in a quasi-random manner: 
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|_:rj designates the fractional part of x, p a and q a are some prime numbers, [Ai,^] and [A' x , A 2 ] are real variational 
intervals which need to be optimized. Parameters Pi and 7$ are obtained in a similar way. The actual values of these 
parameters for the calculation with the largest basis set of 5200 functions is presented in Table I. As is seen from the 
Table fine tuning of variational parameters is not required that greatly facilitates calculations. 

Table II demonstrates a convergence of the variational expansion with the number of basis functions. Extrapolated 
value has been obtained by means of the simple extrapolation formula, 



£(extrap) = E(N) - C x AT 



(3) 



where parameters C and v are taken from the best fit of the last 4 or 5 calculations. 

In Table III a comparison with the most recent and most accurate values is presented. Our result extends accuracy 
of the previous calculations by more than three orders of magnitude. 
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TABLE I: Variational parameters and number of basis functions (rii) for different subsets of the variational wave function with 
N = 5200. Intervals [Ai,^] and [^i,^] correspond to real and imaginary parts of a randomly chosen parameter cti (see 
Eq. (j^) for details), intervals [B\,B 2 ] and [B[,B 2 ] to /3j, and intervals [£71,6/2] and [G'i, G' 2 ] to 74. Prime numbers are p a — 2, 
Pp — 3, p 1 — 5, and q a — 7, qp — 11, q 1 = 13. 





rii 


A 1 


A 2 


A'i 


A' 2 


Bi 


B 2 


Bi 


B' 2 


Gi 


G 2 


G'i 


G' 2 


i = 1 


1160 


0.70 


3.10 


0.00 


0.60 


0.95 


3.35 


0.00 


0.50 


0.00 


0.80 


0.00 


0.75 


i = 2 


1160 


1.00 


9.35 


0.00 


1.05 


0.60 


7.60 


0.00 


0.70 


0.00 


0.70 


0.00 


1.80 


i = 3 


1020 


3.00 


14.0 


0.00 


0.00 


3.50 


14.0 


0.00 


0.00 


0.00 


1.60 


0.00 


1.60 


i = 4 


950 


16.0 


48.0 


0.00 


0.00 


16.0 


48.0 


0.00 


0.00 


0.00 


10.0 


0.00 


5.50 


i = 5 


910 


45.0 


200. 


0.00 


0.00 


45.0 


200. 


0.00 


0.00 


0.00 


30.0 


0.00 


11.0 



TABLE II: Nonrelativistic energies for the ground state of a helium atom °°He. N is the number of basis functions. The last 
digits of the difference in energy between two successive calculations is shown in a third column. 



N 


E nr (in a.u.) 


AE 


3400 


-2.9037243770341195983110931 




3800 


- 2 .903724377034 1195983111421 


490 


4200 


-2.9037243770341195983111540 


119 


4600 


-2.9037243770341195983111572 


32 


5200 


-2.9037243770341195983111587 


15 


extrap 


-2.9037243770341195983111594(4) 





TABLE III: Comparison with other calculations. 



N 



E (a.u.) 



Goldman |7J 

Chuluunbaatar 
et al U 

Sims and 
Hagstrom 

Drake et al. 
This work 



8066 
5669 




-2.903724377034119594 

-2.90372437703411959829 

-2.9037243770341195982999 

-2.903724377034119598305 

-2.903724377034119598311(1) 

-2.903724377034119598311159 



